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This paper demonstrates that a computer aided perturbation theory can easily be realized by use 
of a cumulant approach. In contrast to a recent alternative formulation on the basis of Wegner's 
flow equation method the present approach can be applied to systems with arbitrary Hilbert space. 
In particular an equidistant spectrum of the unperturbed part of the Hamiltonian is not needed. 
The method is illustrated in detail for dimerized and frustrated spin 1/2 chains for which the ground 
state energy is calculated up to seventh order perturbation theory. 
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INTRODUCTION 

Perturbation theory has proved very powerful for the 
investigation of problems which are not exactly solvable. 
Thus, a number of different schemes have been devel- 
oped which are either based on canonical transforma- 
tions (a well-known example is the Schrieffcr- Wolff trans- 
formation 0) or on projection technique Q- Because 
of the increasing computer capacity one would like to 
perform such systematic perturbation expansions up to 
high orders by use of algebraic programmes. Recently 
developed computer aided algorithms are based on the 
flow equation method 0- Q- and Takahashi's formu- 
lation |6| of standard perturbation theory. The latter 
approach was used to derive an effective spin Hamilto- 
nian for high-temperature superconductors Q, whereas 
the flow equation method has been applied to a number 
of low-dimensional systems, see e.g. Refs. 0. However, 
the applicability of the flow equation method is restricted 
to cases for which the unperturbed Hamiltonian has an 
equidistant eigenvalue spectrum. Only then the involved 
set of differential flow equations can be integrated. 

Recently, a systematic perturbation expansion for 
many-particle systems in terms of cumulants has been 
proposed by two of us ||. This method is an projection 
approach and is based on the construction of effective 
Hamiltonians for low-energy properties. Besides the au- 
tomatically preserved size consistency of extensive vari- 
ables, this cumulant method offers compact expressions 
for the different orders of the perturbation theory. The 
main problem for the evaluation of cumulant expressions 
is to properly count the number of contributing configu- 
rations or 'diagrams' to each order. However, since the 
counting of configurations can be taken over by the com- 
puter, the cumulant approach is an ideal starting point 
for an algebraic computer aided evaluation of physical 
quantities. The aim of the present paper is to demon- 
strate how these ideas can be put into practice. 

For the sake of comparability, let us consider a dimer- 



ized and frustrated S = 1/2 spin chain 
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in the limit of strong dimcrization. This model was re- 
cently also investigated in Ref. by use of the flow equa- 
tion method. In Eq. JJJ, (s°) denotes the spin at even 
(odd) site of dimer j. The unperturbed part Ho of the 
Hamiltonian Q describes N uncoupled dimers. There- 
fore, its ground state is the product of singlets on all 
dimers, and the excited states of Ho can be classified with 
respect to the number of local triplets. The perturbation 
Hi describes exchange interactions between neighboring 
dimers. In literature, spin chains of type Q have been 
used for som e sp in-Peierls compounds like CuGeC>3 or 
TTFCuBDT [ljJ[H[Ii- Note that Q is an example for 
a model with an equidistant unperturbed level scheme. 

The paper is organized as follows. In the next sec- 
tion the cumulant approach 9] is formulated. In partic- 
ular, the perturbation expansion of the resulting effective 
Hamiltonian is given in terms of cumulant expressions. In 
Sec. B we develop the computer aided perturbation the- 
ory based on cumulants for the dimerized and frustrated 
5 = 1/2 chain. The ground state energy is calculated up 
to seventh order in the interaction. This model is gener- 
alized in Sec. B so that the spectrum of the unperturbed 
part Ho of the Hamiltonian becomes non-equidistant. Fi- 
nally, the conclusions are presented in Sec. 



CUMULANT APPROACH 

The cumulant approach || starts from the decompo- 
sition of the Hamiltonian H = Ho + Hi into an unper- 
turbed part Ho and into a perturbation Hi . The Hilbert 
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space of the unperturbed Hamiltonian Ho is split into two 
subspaces: The low-energy part Up and the high-energy 
part Uq with projection operators V and Q = 1 — V. 
Let us assume that the two subspaces are separated by 
a finite energy difference. It is our aim to construct an 
effective Hamiltonian for the low-energy subspace Up. 

Motivated by the quantum statistical expression for 
the free energy, the effective Hamiltonian for the subspace 
Up is defined as follows 



H 



off 



ft 



(2) 



where (• • • )p denotes the operator product V{- ■ ■ )V ', ft 
is the inverse temperature Note that due to the 

projectors V in Eq. the effective Hamiltonian H c s 
only acts in the low-energy subspace Up. 

In order to transform Eq. into a cumulant expres- 
sion we introduce generalized cumulants 
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which in contrast to usual cumulants 0, are still 
operator quantities. A detailed discussion of generalized 
cumulants can be found in Ref. 9. By use of series expan- 
sions we transform the effective Hamiltonian (0) into a 
compact cumulant expression. Its Laplace transform can 
be used to derive a perturbation series for the effective 
Hamiltonian. For the case that all states of the relevant 
Up subspace are degenerate with respect to Ho, the re- 
sulting effective Hamiltonian reads at temperature T = 



9) 



H e S (ft -> OO) = 
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Here, Lo is the Liouville operator with respect to Hq. It 
is defined by LqA = [Ho, A] for any operators A. 

To calculate the cumulants in Eq. Q one first decom- 
poses the perturbation Hi into eigenoperators of Lq 



Hi 



with T-iQ%n 



(5) 



Due (J5J also products of % n are eigenoperators of L 
[l^| . Therefore, the energy denominators of Eq. can 
directly be evaluated 
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We are left with the calculation of cumulant expressions 
of the general form (7[ • ■ ■ Tj^)p. For this purpose, we 
rewrite the cumulant expression @ by expanding the 
logarithm into powers of £j and perform the differentia- 
tions. We find the following decomposition of the gener- 
alized cumulants into operator products 
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Thus, the calculation of cumulant expressions is reduced 
to the evaluation of sums over operator products which 
can be easily done by use of a computer. The main limi- 
tation for a concrete realization is given by the increasing 
number of convoluted sums. Consequently, the numeri- 
cal effort may be considerable, if additional restrictions 
have to be taken into account. 



DIMERIZED AND FRUSTRATED SPIN 1/2 
CHAIN 

In this section we show how the cumulant approach 
can be used to perform specific calculations. For this 
purpose we want to construct an effective Hamiltonian 
for the dimerized and frustrated spin 1/2 chain JTJ in 
the limit of strong dimerization. As mentioned above, 
in this limit the unperturbed part Ho of the Hamiltonian 
describes isolated dimers without interaction between dif- 
ferent dimers. Therefore, the low-energy subspace Up is 
given by a single state which is a product state formed 
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TABLE I: Action of the T m (j) as used in the calculations. For 
convenience, the dimer indices of the states are suppressed. 

by singlets on all dimers. The low- and the high-energy 
subspaces are separated by the singlet-triplet splitting 

A = s t - e s = J (8) 

on a single dimer. Thus, the low-energy subspace Up only 
consists of a single state, i.e. the singlet product state. 
The effective Hamiltonian acting in Up can be directly 
identified with the ground-state energy of the complete 

I 
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where Xj-%, 2/3-1, Xj, tjj denote singlets (s) and triplets 
(t~, t°, t + ) on the dimers (j — 1) and j. Let us introduce 
local transition operators %n(j) by the relation 

T-*(j) + r_r(i) + + Ti(j) + = (12) 
= JA {sj • s?_ x + a (s e 3 ■ s]_ x + s° • s^.J } 

which are again classified with respect to the net change 
of the number of triplets. Eqs. (|12|l together with 1)11(1 al- 
lows to construct the transition operators. For instance, 
matrix elements of (jllf) connecting two singulet states 
[on the dimers (j — 1) and j] with two triplet states [on 
the dimers (j — 1) and j] contribute to 7~2(j)- The results 
for T (j), Ti(j), are summarized in Tab. [I] The 

remaining transition operators T_i(j) and T_ 2 (j) follow 
from the relation 71 m (j) = 7^(j). The compact eigen- 
operators T m of the Liouville operator are given by 

N 

Note that the above decomposition of Hi was also derived 
in Ref . H 



problem, multiplied by the projector V. 

Next we decompose the perturbation TL\ into eigenop- 
erators of the Liouville operator Lo As mentioned above, 
the excited states of the unperturbed Hamiltonian 7in 
[see Eq. JQ)] can be classified according to the number 
of local triplets. The creation and the annihilation of a 
local triplet state can be interpreted as the fundamental 
excitation processes. Therefore, also the eigenoperators 
of the Liouville operator can be classified by the number 
of local triplets and TL\ can be rewritten as 

Hi = T- 2 +T_ 1 +%+T 1 +T 2 . (9) 

An eigenoperator T m creates m local triplets. The re- 
spective eigenvalues of the Liouville operator are 

A_ 2 - -2A, A_i - -A, A = 0, (10) 
Ai = A, A 2 = 2A. 

The perturbation Tii [see Eq. JQl] consists of interactions 
between adjacent dimers. Therefore, no further excita- 
tion operators occur, and the T^s can be directly con- 
structed from two-dimer matrix elements 



i 

Now we can calculate the cumulant expressions 

(T mi . . .T mk f p = (14) 

N N 
ji=2 j k =2 

which enter the effective Hamiltonian Q . For that pur- 
pose it is favorable to exploit some additional restrictive 
conditions. In this way the number of contributions to 
the cumulants, which have to be calculated explicitly, can 
be considerably reduced. Remember, the effective Hamil- 
tonian |0J acts in the low-energy subspace Up which is 
formed by the product state of singulets on all dimers. It 
follows 

(i) Products of the T m operators inside the cumulant 
expression do not change the number of triplets, 
i.e. mi + 7722 + • • • + 777fe = has to be fulfilled. 

(ii) Due to the explicit form of the excitation operators 
(see Tab. Pi, the low-energy subspace Up and the 
high-energy subspace Uq are only connected by T 2 
and T- 2 . Thus, non- vanishing contributions to i|14|) 
can only occur if mi = —2 and rrifc = 2. 
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(iii) Due to general properties of cumulants 14] only 
connected processes or 'diagrams' contribute to 
cumulant expressions. Therefore, the maximum 
range of connected dimers in l|14|) is restricted to 



Jmax Jmin 

= max(j 1; 



(15) 



,jk) - min(j'i, . . .,j k ) < {k-1) 



All dimers between j m i n and j n 
lant. 



enter the cumu- 



Whereas (i) and (ii) directly follow from Eqs. I|14|) and 
the condition (iii) represents a basic feature of cu- 
mulant expressions. Note that these restrictive condi- 
tions extremely simplify the evaluation of the effective 
Hamiltonian 0. The zeroth order of H e ff is given by 
7io = Ne s V whereas the first order contribution van- 
ishes due to (Hi)p = 0. The second order contribution 
to the effective Hamiltonian reads 

H e s(/3 -» 00)| 2n dorder = ( 16 ) 

where we have used (ii). Now we insert Eq. 1(1 3() into 
((16(1 and take the condition (iii) into account. Thus 



T~iefi(P — ► OO) 1 2nd order — 
1 N N 

tEE (T- 2 (ji)T 2 (j 2 )) 



(17) 



2A 

]1= 2j 2 =2 



where the dimer site index j can be chosen arbitrarily. 
Note that in 1)17(1 (N — 1) was approximated by N. The 
remaining cumulant in 1(17(1 can be easily evaluated by 
use of Eq. and the matrix elements of Tab. [I] We 
find 



H c ff(/5 — > OO) 1 2nd order 



32 



iVAA 2 (l - 2a) 2 V. 



(18) 



Next, we use the computer to calculate the higher or- 
ders of perturbation theory. By use of the computer 
algebra system Maple V |l6j we have implemented the 
following steps: 

(a) At first, the decomposition JHJ of the perturba- 
tion Tii is inserted in the cumulant expressions so 
that all energy denominators can be easily eval- 
uated [compare with Eq. JJ]J]. Consequently, the 
cumulants of Eq. are decomposed into sums of 
expressions of the form (|14|) . 

(b) Cumulants of the form 1(14(1 are computed by use of 
Eq. if the above conditions (i) - (iii) are fulfilled. 
Otherwise no further evaluation is needed. 



(c) Finally, the remaining operator products of Eq. 
are calculated. 

It is obvious that steps (a) and (b) can easily be imple- 
mented by use of computer loops. For the calculation of 
the operator products of Eq. in step (c), it is sufficient 
to consider finite clusters. The finite order contribution 
of a short-ranged perturbation is independent from the 
cluster size if the cluster is chosen large enough. One can 
also prove explicitly that the results are not affected by 
the cluster size. A state of the cluster is described by 
an array of integers {1, 2, 3, 4} which represent the states 
{s, t + , t°, t~} for each dimer. The T operators applied to 
a cluster change the elements of the integer array which 
then describes the final cluster state. Finally, one has to 
count over all contributing processes. 

As result of the computer aided perturbation theory 
we find the following effective Hamiltonian 



Ti-esi/3 -> oo) 



■-EG)' 



h n V (19) 



The parameters h n up to seventh order read 

, 3 
ho = ^ 

hi = 0, 

h 2 = \{l-2a)\ 

h 3 = + U-2a) 2 , 



13 27 3 2 \ M n , 2 



(20) 



+ Y a--a' ) (l-2a) z , 



, 89 311 93 2 _ o \ , „ n2 

h*, = 1 a H a 2 - 45a 3 (1 - 2a) , 

24 12 2 1 



/463 454 1307 2 

1 n A n, z — 



a A - 84a J - 159a 4 x 



h 7 = 



V 48 9 6 

x (l-2a) 2 , 
/81557 257909 215995 , 173579 



V 3456 1728 



432 



216 



14865 4 879 B Wi „ , 2 
+ — «° ) (1 - 2a) 2 



Apart from the projector V the effective Hamiltonian 
[ (imp and 1(20(1] can be directly identified as the ground- 
state energy of the original Hamiltonian . As discussed 
above, this follows from the fact that the low-energy sub- 
space Up only consists of a single state. Note that the 
present result for the ground-state energy agrees with the 
result obtained before on the basis of the flow equation 
method 0. 
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THE GENERALIZED MODEL 

The unperturbed part Ho of the Hamiltonian ijljl. 
which we have studied up to now, has had an equidis- 
tant spectrum. Note that this property of Ho was ex- 
plicitly needed in the flow-equation approach to construct 
perturbative effective Hamiltonians |8j- In the following 
we want to demonstrate that this property of the model 
Hamiltonian is not required in the cumulant method dis- 
cussed above. For that purpose let us modify the unper- 
turbed Hamiltonian Ho of the model JD, whereas Hi is 
kept unchanged: The coupling between the two spins of 
each dimcr is chosen to be different for dimers with even 
and odd site index j 

N N/2 

Ho = J s r s ° + JS J2 s h ■ s h • ( 21 ) 

4=1 4=1 
Note that the new Hamiltonian Ho does not change the 
eigenstates of the previous unperturbed Hamiltonian Q . 
However, the eigenenergies are different. The dimer sin- 
glet and triplet energies depend on the dimer index j and 
are given by 

e s j odd 

(1 + 8)e s j even 



sti) = 



(22) 



(23) 



e t j odd 

(1 + 6) et j even 

Therefore, the spectrum of the unperturbed Hamiltonian 
Ho is not equidistant anymore. Due to the modification 
(|2*T|l of Ho also the Liouville operator L has changed. 
The decomposition of the perturbation Hi, Eq. (JJJ, into 
transition operators with respect to Ho has to be mod- 
ified too. In particular, the number of created local 
triplets can not be the only classification criterion any- 
more since the singlet-triplet splittings differs for different 
dimer sites. In fact, the former transition operators Ti, 
Tq, and 71 1 have to be split up. The perturbation Hi 
can now be written as 

Hi = T-2 + 71i,_ + 71 1)0 + To.- + To.o (24) 
+ To.+ + Ti fi + Ti t+ + T 2 

and the respective eigenvalues of the modified Liouville 
operator now read 

(25) 



A 2 


= {2 + S)A 


Ai,+ 


= (1 + 5) A 


Ai.o 


= A 


A ,+ 


= 6A 


A ,o 


= 


Ao,- 


= -6 A 


A_i,o 


= -A 


A-i,- 


= -(1 + 6)A 


A_ 2 


= -(2 + 6) A 



The second index x in T±i^ x and T), x denotes how 
the number of triplets on dimers with even index j is 
changed, (x — +(—) describes the creation (annihila- 
tion) of a triplet and x = no change). Note that the 
original model Q is given by 6 = 0. 

The next steps can be done as before. At first we in- 
troduce local transition operators. Furthermore, the re- 
strictive conditions (a), (b), (c) are still valid so that the 
evaluation of the cumulants can be done similar as be- 
fore. (Of course, now one has to distinguish between odd 
and even dimer indices.) By use of the computer aided 
perturbation theory we now find an effective Hamiltonian 
which has the form of expression (|lfJ|) . The parameters 
h n , up to fourth order, now read 



h = 

hi - 

h 2 = 

h 3 = 

hi = 



;(2 + <5), 



(26) 



0, 



1 (1 - 2a) 2 , 



(2 + 6) 
4 

(2 + 5) 



2 y -+3a)(l-2aY 



3^^112 



2 + 6 



(2 + S) 2 V (3 + 2<5)(3 + <5) 

(2 + 6){l- 26) 2 5 - 52a + 20a 2 



1 + 6 



2 + 6 



In fact, we have calculated the h n also up to seventh 
order. Since the expressions are rather involved, here the 
h n 's are only given up to the fourth order. Higher orders 
are available on request. Note that in the case of 6 = 
Eq. reduces to (|5Ujl . 



CONCLUSIONS 

In this paper we have shown that a recently developed 
cumulant method can be used to develop computer 
aided perturbation theory. Size consistency of extensive 
variables is fulfilled. The cumulant method offers com- 
pact expressions for the different orders of the perturba- 
tion theory. Furthermore, the evaluation of the cumulant 
expressions is reduced to the problem how to properly 
count the contributing processes. 

We have applied the cumulant method to the dimcr- 
ized and frustrated spin 1/2 chain. For this model the 
ground-state energy was calculated up to seventh order 
perturbation theory where our results agree with those 
obtained by the flow equation method |8|]. It turned out 
that the efficiency of the computer based evaluation of 
the cumulant expressions can be enormously improved 
if restrictive conditions are considered. In this way the 
number of vanishing contributions in the calculations can 
be reduced. Furthermore, we have modified the model by 
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an additional site-oscillating dimer coupling so that the 
spectrum of the unperturbed part of the Hamiltonian is 
not equidistant anymore. For the generalized model we 
have calculated the ground-state energy up to seventh 
order perturbation theory as well. Note that in contrast 
to the cumulant approach the flow equation method used 
in Ref. 8 requires an equidistant spectrum of the unper- 
turbed part of the Hamiltonian. On the other hand, the 
flow equation method can also be used to calculate exci- 
tation energies. 

The derivation of the cumulant expression jlj for the 
effective Hamiltonian H e g(/3 —> oo) was based on the 
assumption that it is acting in the lowest energy sub- 
space Up of Ho ■ This subspace can either be degenerate 
or nondegenerate. For the dimerized spin 1/2 chain the 
lowest energy subspace is one-dimensional and is given 
by the product of all dimer singlet states. An example 
for a degenerate unperturbed groundstate was discussed 
in Ref. 9 1 . Note that the cumulant expression (@J is also 
closely related to an effective Hamiltonian which was de- 
rived by Takahashi Q. There, the eigenvalue problem 
of the full Hamiltonian is transformed to that of an ef- 
fective Hamiltonian which acts in a degenerate or non- 
degenerate energy subspace of Hq. By comparing the 
perturbation expansion order by order with that of 
Ref. 6] the equivalence of both approaches can be shown. 
Takahashi's approach does not involve the temperature. 
Therefore, the subspace in which the effective Hamilto- 
nian acts is not necessarily the lowest energy subspace 
of Hq. Thus, one might expect the cumulant result (@J 
should also be valid for this case which would allow to 
calculate excitation energies. Finally, a projector-based 
renormalization method (PRM) for effective Hamiltoni- 
ans was recently introduced by two of the present authors 
Fjj. By using perturbation theory also in this approach 
18j a close relation to the cumulant expression Q can be 
found. However, the PRM treatment seems to be more 
suited to calculate excitation energies than the cumulant 
approach presented here. 
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